Proximity-induced superconductivity in graphene 
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We propose a way of making graphene superconductive by putting on it small superconductive 
islands which cover a tiny fraction of graphene area. We show that the critical temperature, T c , can 
reach several Kelvins at the experimentally accessible range of parameters. At low temperatures, 
T <C T c , and zero magnetic field, the density of states is characterized by a small gap E g < T c 
resulting from the collective proximity effect. Transverse magnetic field H g (T) oc E g is expected to 
destroy the spectral gap driving graphene layer to a kind of a superconductive glass state. Melting 
of the glass state into a metal occurs at a higher field H g 2(T). 

PACS numbers: 74.78.-w, 74.20.-z, 74.81.-g 



Among numerous fascinating properties, graphene [l], 
2] provides a unique possibility to study the phenomenon 
of proximity-induced superconductivity in very favorable 
conditions. Experimental studies of the Josephson cur- 
rent through graphene in standard wide planar SNS junc- 
tions Q have shown that proximity effect in graphene is 
qualitatively similar to the one known for usual dirty 
metals. In this Letter, we show that even a tiny amount 
of graphene area covered by small superconductive is- 
lands (with good electric contact to graphene) can lead to 
a macroscopically superconductive state of the graphene 
film, with T c in the Kelvin range. 

We consider a system of superconductive (SC) islands 
of radius a (with the typical value of few tenths of 
nanometer) placed approximately uniformly on top of a 
graphene layer (with the typical distance between islands 
b in the sub-micron range) shown in Fig. [TJ We assume 
that b is much larger than both a and the graphene mean- 
free-path I. Moreover, present theory will be limited by 
the case I ^ a when electron motion in graphene is dif- 
fusive at all relevant scales. We will not be particulary 
interested in phenomena in the vicinity of the graphene 
neutral point, assuming relatively large gate potentials 
\V g \ > 10 V, and carrier density n > 10 12 cm 2 . We as- 
sume graphene Fermi energy Ep 3> Ao 3> T c , where Ao 
is the island's superconductive gap. Graphene sheet can 
be either single- or few-layered: the only relevant features 
are (i) high diffusion constant D > 10 2 cm 2 /s, and (ii) 
very low (in comparison with metals) electron density, 
which allows to combine moderate values of dimension- 
less conductance g — (h/e 2 Ru) > 3 with high Thouless 
energy i?Th = hD/b 2 . Not very large values of sheet 
conductance g are practically favorable to avoid suppres- 
sion of superconductivity in small SC islands due to the 
inverse proximity effect 

Below we treat graphene as a normal diffusive 2D metal 
within the standard approach based on the Usadel equa- 
tion [|| ; its applicability to diffusive graphene was proven 
m Ref. 0,0]. The intrinsic Cooper channel interaction 
in graphene can be neglected due to its low DOS 



FIG. 1: (Color online) Graphene film covered by supercon- 
ducting islands. 



Similarly, phonon-induced attraction is also weak. 

Proximity coupling and transition temperature. — We 
start with calculating the Josephson coupling energy be- 
tween two superconductive islands of radius a separated 
by distance 6 » o, neglecting the presence of other is- 
lands. Such a pair-wise approximation is adequate for 
determination of T c , but breaks down at T < T c j ln(6/a), 
as shown below. For the single SC island on graphene, 
the Matsubara-space Usadel equation for the spectral an- 
gle 9 U and corresponding boundary conditions read as 



L>V 2 (L - 2Msin<L = 0, 



dr 



27ra 



■ cos 6 U 



= 0. 
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The normal (G u ) and anomalous (F u ) components of 
the matrix Green function the in Nambu-Gorkov space 
are expressed via the spectral angle 8^ and the order 
parameter phase <p as G„(r) = cos6> w (r) and F u (r) = 
e lip sin B u (r) . The full matrix structure of the anoma- 
lous Green function with the valley and spin spaces 
included (Pauli matrices tt and s, respectively) is deter- 
mined by the usual s-wave pairing in the SC islands: 
Fui oc Tx x s y . The interface conductance G- lnt is treated 
below as a phenomenological parameter which accounts 
for Fermi velocity mismatch and a potential barrier on 
graphene- metal interface [l0( . 

It is crucial for further analysis that the two- island gen- 
eralization of the nonlinear problem ([T]), (|2|) can be lin- 
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earized while calculating the Josephson current at inter- 
island distances b ^> a. Indeed, the total current can 
be calculated by integrating the current density over 
the middle line between the islands, on the distance 
Pi, 2 > b/2 from them. This procedure also involves 
summation over Matsubara energies u n = nT(2n + 1), 
with the major contribution to the sum coming from 
w n ~ i?Th- At such uj n and the spectral angle 9 
is small, and linearization of Eqs. |T]) and |2|) leads to 
the solution 



,(r)=A(cj)K 



0(t u 



hi(L u /a) ' 



(3) 



with = y/D/2u, t w = (G int /27r 5 ) In(L w /o), and 0(i) 
solving the equation Q(t) = t cos O(t). The function 
A(ui) evolves between the tunnel and diffusive limits as 



GW(27T fl ), 

w/[2ln{L u /d)], 



G int < 2irgfhx{L u /d), 
G int > 2wg/\n(L u /d), 



(4) 



and is always small for ln(L w /a) ^> 1. Thus the Joseph- 
son current I(ip) = I c sin ip between two SC islands with 
different phases, ipi~(f2 — <p, can be calculated using the 
linearized two-island solution for the anomalous Green 
function: F u (r) = e^ 1 sin0 w (|r - n|) + e^ 2 sin0 w (|r - 
r 2 |). The standard calculation of the Josephson energy 
Ej = (h/2e)I c then leads to 



Ej(b, T) = AirgT £ A 2 (uj n ) P(y/u n /83r h ) 



(5) 



w„>0 



where P(z) = z Kq(z cosh t)K\{z cosh t) dt. 

A two-dimensional array of SC islands with the cou- 
pling energies ([5]) undergoes the Berezinsky-Kosterlitz- 
Thouless transition at 



T c = jEj(b,T c ), 



(6) 



where the numerical coefficient 7 depends on the array 
structure. Below we will assume that the SC islands 
form a triangular lattice, in which case 7 ~ 1.47 11|. 
For the interface conductance G; nt comparable with the 
sheet conductance g, one finds the transition temperature 
T c ~ £rh- I n general, T c can be obtained by numerical 
solution of Eq. © using Eqs. © and (0. The result ob- 
tained for the ratio T c j E^ as a function of Gi n t for g = 6 
(R n ~ 700 0) and 6/a = 10 is presented in Fig.H With 
the graphene diffusion constant = 500 cm 2 /s (see, e.g., 
(l^ |) and & = 0.5 jum, one estimates i?xh ~ 1-5X, leading 
to T c in the range 1 3 K for 5 < G in t < 20. 

Low temperatures: spectral gap and order parameter. 
Now we switch to the low-temperature range T -C T c and 
consider the issue of the spectral gap for the excitation 
above the fully coherent ground state (with all phases ipi 
equal). The density of states v{E) — vq Re cos 9(E) is de- 
termined then by the periodic solution of Eqs. ([T]) and 
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FIG. 2: The critical temperature, T c , and the zero- 
temperature spectral gap, E g , vs. the interface conductance 
Gint (the sheet conductance g — 6, and b/a — 10). 



analytically continued to real energies: \ui\ — > iE. This 
periodic problem is equivalent to the one defined within 
the single (hexagonal) elementary cell, supplemented by 



the additional condition riV9\ 



rer 



0, where T is the cell 



boundary. Solution of the Usadel equations for such a ge- 
ometry leads to formation of the spectral gap E g similar 
to the minigap for one-dimensional SNS junctions [l3j ]. 
To find it, we write 0(r) = 7r/2 + iip(r) and determine 
the spectral boundary as the value of E where equation 



2£cosh^ = 



(7) 



ceases to have solutions with real ip(r) [14|. At large 
ln(6/a) one may approximate the hexagonal boundary T 
of the elementary cell by the circle of radius R = b/2. 
For the ideally transparent interface, numerical solution 
of the radially symmetric equation ([7]) gives for the value 
of the zero-temperature spectral gap: 



E„ 



HD/R 2 



2.63 



Th 



1.52 ln(R/a) - 1.2 m(&/4a)' 



(8) 



Decreasing the interface conductance Gi„t leads to the 
suppression of the the minigap, as shown in Fig. [2j 

In the limit of large ln(6/a), the spectral gap E g <§; T c . 
Smallness of the gap distinguishes the system with su- 
perconductive islands from usual dirty superconductors. 
Roughly speaking, it behaves as a continuous 2D super- 
conductor at the energy/temperature scales smaller than 
Eg, whereas in the range E g < (E, T) < T c it can rather 
be described as an array of weak Josephson junctions. 

The existence of the sharp gap ([8|) in the electron spec- 
trum looks surprising, as only a tiny fraction (a/6) 2 of 
graphene area is in direct contact with SC islands. The 
presence of this gap can be traced back to the periodic 
structure of islands we assumed. Therefore any irregular- 
ity in the positions of SC islands will lead to the smear- 
ing of the hard gap. Assuming that islands' locations are 
shifted at random from the sites of the ideal triangular 
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lattice, with the typical shift 6b <C b, one can reduce the 
problem to the effective one, defined on a scales large 
than array lattice constant. Random displacements of 
islands will be seen, in terms of this effective model, as 
local fluctuations of the superconductive coupling con- 
stant HI, lj| , leading to the smearing of the gap with the 
relative width 8E g ~ (5b/b) 2 E g . The sharp gap will also 
be smeared by thermal fluctuations of island's phases tpi 
and finite thermal coherence length Lt- Thus we expect 
the spectral gap to be observable at T <C E g . 

Even in the presence of the gap smearing, strong sup- 
pression of the local DoS in graphene at E < E g , should 
be seen by the low-temperature Scanning Tunnelling Mi- 
croscopy. The spectral (pseudo) gap is a signature of col- 
lective proximity effect which cannot be quantitatively 
described by a pair-wise interaction between SC islands 
as soon as low-energy scales < E g are involved. The 
corresponding spatial scale 



ig = \ KD/E g ~ b^Hb/a) 



(9) 



plays the role of low-temperature coherence length in the 
(dirty-limit) superconductor. Under our main condition 
ln(6/a) ^> 1, the coherence length £ g 3> b, which allows 
continuous treatment of the array at low temperatures. 

The local superconductive order parameter in 
graphene, T{v) = J dujF w (r), can be found at T < E g as 



(r-rjflnQr-Tjl/ay 



(10) 



where Yj are the coordinates of SC islands, t{r) = 
(Gi n t/2ng)ln(r/a), and we used the solution © for 
u) > E g . The divergent sum in Eq. (JTUJ) should be cut 
at |r — Tj \ ~ £ g since the spectral gap © suppresses the 
lowest-w contribution to ^(r). Equation pop is not ap- 
plicable in the vicinity of SC islands since Eqs. © and 
© cannot be linearized at small |r — r_j|. 

At zero magnetic field, ipj — const and the space- 
averaged order parameter T is given by 



T 



d 2 r De\t(rl 
r 2 ln(r/a) 



ir 2 lnln(&/a) , , 
= TT-riiD - ,) ' (11) 
2 \n(b/a) ' v ; 



where the last expression refers to the large-Gi n t limit 
(O ps 7r/2) and n; i=s 1/b 2 is the concentration of SC is- 
lands. Comparison of (TTTj) and © provides the condition 
for neglecting the intrinsic Cooper-channel interaction in 
graphene, \ g . Namely, its presence would generate the 
the energy gap A g = \ g T. This "intrinsic" gap can be 
neglected compared to proximity-induced gap ([5]) pro- 
vided that A g <C 0.5. Comparison with the estimate for 
intrinsic Cooper interaction constant |8y shows that the 
latter is indeed negligible. 

Electromagnetic response. — Linear response of the su- 
perconductive film to a weak electromagnetic field is 



characterized by the superconductive density p s . In the 
intermediate temperature range, E g <C T <C T c , one can 
easily calculate p s within the pair-wise approximation for 
the proximity coupling: 



p. CO 



-f 27rr 3 Ej(r,T)dr. (12) 

Taking Ej(r,T) from Eq. © wc find 

p,(T) = (7T 3 /3)gA 2 (nT)E 2 jT, (13) 

where g = 2vqD and the numerical factor corresponds 
to the triangular array with m = (2/V3)b- 2 . At lowest 
temperatures, T < E g , the function p s {T) saturates at 
the value that can be estimated by the replacement T — > 
E g in Eq. (|13[) . For highly transparent interface we obtain 



p.(0) ps 10 ffj B T h/ln(6/a). 



(14) 



Comparing (p~4|) and (jSJ) we find that p s (0) ~ ^gEg, which 
is typical for dirty superconductors with the gap E g . 

The critical current density per unit length at the low- 
est temperatures, T <C E g , can be estimated as 



2e 7r^ cqD 



(15) 



where we used, as an estimate, the T = limit of the 
pair- wise Josephson coupling energy © at G; nt /g ^> 1: 



Ej(b,0) 



gE 



Th 



In (b/a) 



(16) 



The effect of the transverse magnetic field is character- 
ized by two different field scales: 



H„ 



2< 2 



0.4 $ 
ln(6/4a) ~W'' 



H, 



Glass 



b 2 ' 



(17) 



In the low-field, low-T region (B <C H g and T < E g ), 
magnetic field produces well-separated pancake "hyper- 
vortices" with the core size £ g > b (the local DoS is 
gapless in the core regions). These vortices are strongly 
pinned by the underlying array structure, so a high crit- 
ical current j c \ ~ jc(0)(fe/£s) is expected. At B w H g 
vortex cores overlap and the proximity gap is totally de- 
stroyed, so Hg is an analogue of the upper critical field 
H C 2- However, the metallic state is not formed right 
above H g , at least at T -C E g . In this field range 
one deals with a system of frustrated pair-wise Joseph- 
son couplings, with full frustration achieved at B » 
Hciass- In this high- field range, average values of Joseph- 
son coupling are exponentially suppressed, Ej(B) oc 
exp(— B/Hqiuss). However, as shown in [l?], EH, actual 
(random-sign) Josephson couplings are much stronger 
due to mesoscopic fluctuations: 



\b) = [{Ej(b)f 



1/2 



ln 2 (6/a)' 



(18) 
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H I \ 

strongly 
frustrated 
SC glass 

H G 1 * 



weakly frustrated SC glass 



inhomogeneous ' 
bulkSC 



Josephson-coupled islands 



E a 



T c T 



lie state [19j. The same problem of suppressed Andreev 
conductance appears to be even more serious with the 
decrease of the electron density towards the graphene 
neutral point: scattering cross-section of electrons on the 
SC islands drops in the range of kpa ~ 1 leading to effec- 
tive decoupling of island's phases and to strong quantum 
fluctuations. These fluctuations ma y le ad to the QPT of 
the superconductor-metal type 1|| 2oj |. 

We are grateful to P. A. Ioselevich, P. M. Ostrovsky 
and M. Titov for useful discussions. This work was par- 
tially supported by RFBR Grant No. 07-02-00310. 



FIG. 3: Schematic phase diagram of the graphene sheet with 
superconductive islands; all lines refer to crossovers rather 
than to sharp phase transitions. 



which is just by the factor 1 jg smaller than the (T, B) = 
pairwise coupling (|16[) . The estimate (fT5|) shows that at 
T — > the superconductive glass state survives up to high 
magnetic fields H g2 (0) S> -ffciass • The value^of H„^ (0) is 



determined by quantum phase fluctuations [18 



H, 



92 



2 VA 



(19) 



The overall phase diagram in the (H, T) plane is shown 
schematically in Fig. [31 We emphasize that the lines in- 
dicated do not refer to sharp phase transitions (which are 
absent in the presence of magnetic field, apart from some 
special values of frustration for the case of well-defined 
lattice of islands) but rather mark a crossover regions. 
Note that determination of T c for a rational values of frus- 
tration / = |, |, § (where results for nearest- neighbor 
XY model are available) is complicated by the necessity 
of accounting for long-range proximity couplings. 

To conclude, we have shown that graphene can be 
made superconductive with T c of the order of few Kelvins, 
due to collective proximity effect induced by small super- 
conductive islands covering only a tiny part of graphene 
sheet area. The spectral gap is expected at low temper- 
atures T <C T c and magnetic fields B <C $o/b 2 . Trans- 
formation from a continuous disordered superconductive 
state to a weakly coupled junction array is predicted with 
the temperature and/or magnetic field increase. 

Our study was based on the standard Usadel equa- 
tions which are valid for sufficiently disordered samples 
(I « a). In the opposite limit of quasi-ballistic electron 
motion around islands, the Andreev subgap conductance 
might decrease, leading to the increase of quantum fluc- 
tuations of phases tpi (t) [l|J [2(| • Quantum fluctuations 
can be neglected under the condition b 2 \n(b/a) <C b 2 cl 
where b c is the critical distance between the islands mark- 
ing the quantum phase transition (QPT) to the metal- 
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